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The constraints imposed by nano- and microscale confinement on the conformational degrees of 
freedom of thermally fluctuating biopolymers are utilized in contemporary nano-devices to specifi- 
cally elongate and manipulate single chains. A thorough theoretical understanding and quantifica- 
tion of the statistical conformations of confined polymer chains is thus a central concern in polymer 
physics. We present an analytical calculation of the radial distribution function of harmonically con- 
fined semiflexible polymers in the weakly bending limit. Special emphasis has been put on a proper 
treatment of global modes, i. e. the possibility of the chain to perform global movements within the 
channel. We show that the effect of these global modes significantly impacts the chain statistics 
in cases of weak and intermediate confinement. Comparing our analytical model to numerical data 
from Monte Carlo simulations we find excellent agreement over a broad range of parameters. 



PACS numbers: 87.15.ad,82.35.Lr,36.20.Ey 

Understanding the statistical properties of biopoly- 
mers on the single molecule level has been one ma- 
jor objective of polymer research over the past decades. 
While our theoretical knowledge about semiflexible chain 
molecules in free space is developed to a fairly sophis- 
ticated level, a steadily growing theoretical interest in 
the mechanical properties of semiflexible polymer solu- 
tions [1-5J, and networks [SHU], as well as the impe- 
tus from nano-technological applications [T^l [T3J calls for 
an expansion of this knowledge to the case of polymers 
trapped within some confining environment. Odijk |14j 
and de Gennes |15j laid the groundwork for further an- 
alytical treatments of confined chain molecules by pro- 
viding very successful scaling theories for confined semi- 
flexible (Odijk) and flexible (de Gennes) chains. Recent 
advances in photolithographic techniques allowed for the 
fabrication of channel-like structures down to nanome- 
ter scales, thus enabling the experimentalist to investi- 
gate the effect of confinement on single chain molecules. 
Channels of different sizes have been used to directly ob- 
serve the fluctuations of F-actin and DNA molecules and 
to investigate the respective statistical properties within 
Odijk's 16 18J and de Gennes' [13] scaling regime and 
at the crossover between both [T2] . 

On the theoretical side, questions related to thermo- 
dynamic properties and the end-to-end distance statistics 
of polymers in confining geometries have largely been ad- 
dressed by means of computer simulations [TM2"!2] . A re- 
cent contribution to analytically capture the radial dis- 
tribution function (RDF) [23] of semiflexible chains in 
cylindrical confinements was based on the physical idea 
that the polymer's transverse fluctuations are effectively 
suppressed along the contour, if the confining channel 
is sufficiently narrow. In order to facilitate calculations, 
this approach then employed torqued ends boundary con- 
ditions [23] within a channel- fixed reference frame [25]. 



assuming that the particular choice of boundary condi- 
tions is of minor importance for the determination of the 
chain's end-to-end distance statistics. While feasible in 
a strong confinement regime, where the polymer almost 
perfectly aligns with the symmetry axis of the channel, 
the chain statistics predicted by this approach is at vari- 
ance with our results from Monte Carlo simulations for 
intermediate and weak confinement strengths. 

In this letter we go beyond this strong confinement ap- 
proximation and provide a thorough analysis of the end- 
to-end distance statistics of semiflexible chains, trapped 
in cylindrically symmetric confinements, based on ana- 
lytical and numerical methods. Special attention is paid 
to the role of global chain movements, which is techni- 
cally captured by means of a careful choice of boundary 
conditions. As a central result we will present a refined 
version of the RDF which is shown to reproduce the chain 
statistics even in cases of weak and intermediate confine- 
ment. In contrast to the RDF obtained in ref. [53], our 
approach is capable of reproducing the correct form of the 
RDF in the particular limit of zero confinement [JSJ [26] . 
In order to test our analytical predictions we employed a 
standard Monte Carlo (MC) scheme [T5]. For this pur- 
pose we used a discretized version of the wormlike chain 
model, consisting of N inextensible segments of length 
N^ 1 (where the contour length L = 1 is set to unity). An 
energetic penalty, proportional to the squared transverse 
displacement from the symmetry axis of the cylindrical 
potential, was calculated at the ends of each segment 
to mimic the presence of confining channel walls. Self- 
avoidance effects, which are not important for sufficiently 
stiff chains, have been neglected. During the simulation 
process the chain's ends were assumed to be completely 
free. 

To reveal the importance of boundary effects in the 
context of end-to-end distance statistics, we start with a 
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FIG. 1: Definition of polymer- and channel fixed reference 
frame £ p (red, dashed) and E c (black, solid). 9 measures the 
angle between the two z-axes. 



ter c = L^ 1 may be pictured as the average number of 
collisions of the contour with the channel walls. 

In order to proceed analytically, we investigate eq. ([IJ) 
in Laplace space, i.e. we calculate a generating function 

£(0 = £ [G{r)\ ex / V[v{s)\ exp (-£[r(s), £]) , (4) 



where £ [.] = / °° dpe~^ p (.) denotes the Laplace trans- 
form with respect to the stored length p = 1 — r. Thence 
the RDF itself is given by 



short review of some qualitative features concerning the 
connection between the chain's shape and the constraints 
at its boundaries. Reference [T5] gives a classification of 
the shapes, semiflexible polymers in intermediate con- 
finements actually attain. There it was shown that the 
distribution of transverse displacements for chains whose 
end-monomers are free to move within the channel resem- 
bles a bone-like shape with particularly enhanced fluctu- 
ations within a boundary layer of size ~ Ld (L^'. Odijk's 
deflection length |14j). In contrast, hinged ends bound- 
ary conditions give rise to qualitatively different, cigar- 
like shapes of the polymer. This, as will be detailed in 
the course of our quantitative discussions, results in sig- 
nificantly larger end-to-end distances in the intermediate 
confinement regime and significantly smaller end-to-end 
distances in the limit of vanishing confinement respec- 
tively. In the limit of strong confinement the chain's 
shape approaches a straight line in both cases and the 
quantitative differences related to the end-to-end dis- 
tance statistics vanish. 

The RDF of the system may be determined by means 
of an integral over all configurations r(s) consistent with 
one particular value r of the end-to-end distance R, and 
weighted by the appropriate Boltzmann factor: 

G(r) oc J V[r(s)}6 (R[r(s)} - r) exp {-U[r{s)]) . (1) 

Here and in the following all energies are measured in 
units of ksT and all lengths are measured in units of 
the contour length L. The system's Hamiltonian H = 
Wbnd + "Hpot consists of two parts. The bending part is 
given by a standard wormlike chain Hamiltonian \27\ in 
its continuous form HE\ 
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where e, the inverse persistence length, measures the 
chain's flexibility. The second part accounts for the chan- 
nel confinement, which is modeled by a harmonic poten- 
tial 
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Here r_L denotes the transverse displacements from the 
symmetry axis of the channel and the collision parame- 



G(r) = [£(0] = 



d^ 
2ni 



(5) 



In the limit of strong confinements, where the polymer's 
transverse fluctuations from the channel axis become suf- 
ficiently small (|r'jj 2 = 0(e/c) <C 1), one may track the 
polymer's position with respect to the channel axis and 
employ a weakly bending rod ( WBR) approximation |25j . 
Then, the exponential £ attains the following form 



£[r( S ),£]=~{rlr'j!, 



(6) 



1/2 

where (x|y) = J -^^ds [x'(s)y(s)] abbreviates the L2 

scalar product ("x*" meaning "transpose of x") and 
where primes indicate derivatives with respect to s. Note, 
in particular, that within the WBR approximation the 
exponential £ becomes a functional of the transverse dis- 
placements r±(s) only. As shortly discussed above, the 
particular choice of boundary conditions is of minor rel- 
evance in that limit. One may, therefore, make the sim- 
plest choice of either hinged or torqued ends boundary 
conditions in order to make the surface terms in the first 
line of eq. ([6| vanish. This then allows to use sim- 
ple trigonometric eigenfunctions which diagonalize the 
operator occurring in the second line. Both choices fi- 
nally lead to the strong confinement approximation of 
the RDF, determined in ref. |23j . 

It must be stressed, however, that the particular form 
of the exponential £ stated in eq. ^ certainly brakes 
down in the limit of vanishing confinement, where the 
transverse fluctuations rj^ from the channel axis may be- 
come large and the WBR approximation fails. Moreover, 
as detailed above, boundary effects become important 
outside the strong confinement regime, whence the some- 
what artificial choice of hinged or torqued ends boundary 
conditions within a channel-fixed reference frame is no 
longer appropriate. Both problems can be remedied by 
a proper choice of coordinates. To this end, in addi- 
tion to the channel-fixed reference frame S c used so far 
(whose z-axis coincides with the symmetry axis of the 
system), we define a second "polymer- fixed" reference 
frame E p , whose z p -axis shall be defined to connect both 
end-monomers of the polymer (cf. fig. |T|) [30 . Within 
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this newly denned frame we may apply Monge parame- 
terization 



r P 0) = 



?±A S ) 
s + 0{a) 



s e 



1 1 

2' 2 



(7) 



where a — e/c for intermediate (and strong) confinement 
(c > 1) and a = e for vanishing confinement (c — > 0). 
Note that the physically correct boundary conditions 



within E p read 



r -L,p(»| 



= ±1/2 
= ±1/2 



= 



(8a) 
(8b) 



Here eq. ( 8a ) is a consequence of the very definition of 



the polymer-fixed reference frame, whereas ( 8b I arises 



from the physical requirement of a vanishing torque at 
both ends of the polymer. Within this formulation of 
the problem the WBR approximation holds true when- 
ever a < 1, which may be realized either by virtue of 
sufficiently strong confinement strengths c or by virtue 
of sufficiently small chain flexibilities e. The exponential 
in eq. Q then takes the form 

£[r p (s),£,0,u] = 

2 e | r " ,p r ±,pl-i — r T,p r ±,p|l i } + 2 r J-.p rj ->p|- 



(9) 



where Pj_r p = rj_,p and 
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cos 9 sin ( 
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(10) 



gives the transverse displacement vector with respect to 
the channel- fixed frame E c . Here 9 is the angle between 
the z-axes of the two reference frames, and u = (u x , 
connects the origins of E c 
the transformation in eq. 
transverse polymer fluctuations rj_ p in terms of simple 
sine-modes 



and E p . Figure [T] illustrates 



(10). An expansion of the 
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renders the exponential in eq. (|9| Gaussian and, in ad- 
dition, is consistent with the boundary conditions within 
E p , as stated in eqs. ([8|. 

Note, that we were able to preserve hinged ends bound- 
ary conditions and all the convenient features brought 
about by these at the expense of dealing with two refer- 
ence frames, rather than a single one. Technically this 
means, that we have to deal with additional degrees of 
freedom, embodied by the rotation angle 9 and the dis- 
placement vector u, which would be absent if we stipu- 
lated hinged or torqued ends boundary conditions within 
E c . There is a conceptual difference between the lat- 
ter degrees of freedom and the degrees of freedom rep- 
resented by the expansion coefficients {a„} in eq. (111. 



While variation of the expansion coefficients {a n } at con- 
stant values of (u, 9) leads to deformations of the contour 
at fixed relative position of the two sets of axes E c and 
E p , altering the values of (u, 9) while fixing {a„} amounts 
to shifting a frozen contour with respect to the channel. 
Therefore, in what follows we will refer to the coefficients 
{a n } as "internal modes" and we will reserve the term 
"global modes" for the entities (u, 9). In the following, 
this conceptual difference will help us illuminate the dif- 
ferent statistics due to different boundary conditions. 

For the time being we concentrate on the internal 
modes, treating the global modes as parameters. We thus 
calculate 
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Z 



exp (-£[{a„},£,0,u]) , 



(12) 

which is the Laplace transform of the joint probability 
distribution function Q{r, 9,u). The integrals occurring 
in eq. ( 12 ) are Gaussian and may readily be performed. 



It should be noted, however, that the resulting expression 
is a rather involved function of the global modes, which is 
why we resort to the following expansion scheme in order 
to render upcoming cos 9- and u-integrations amenable 
to analytical calculations: 



C(^9,u) = ^e 



- n2 * e -^ u2 LW(0x 
l + sin 2 + L«(e)u 2 



(13) 



Here we have chosen to perform an expansion with re- 
spect to the small quantities sin 2 9 and u 2 in order to 
account for the symmetry properties of the system un- 
der study. Also note that the expansion in eq. (13) 
explicitly keeps track of the rigid rod limit e —¥ 0, where 

C oc e~5K sln2 e e~fc u2 . In what follows we will calculate 
the effect of the global modes to leading order. 

Using eqs. ^ and (12), the zeroth and first order 
coefficients in expansion ( 13 ) are found to be 
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Comparison of Q{r) (eq. \2Q\ . solid lines) and Qh{r) 
dashed lines) with MC simulation data (data 
points) for weak and intermediate confinement strengths, (a) 
Weak confinement (c = 0.5, e = 0.7). For entropic reasons the 
global modes shift the end-to-end distance statistics towards 
larger values of r. (b) Intermediate confinement (c = 3.0, 
e — 0.5). The energy optimizing effect of the global modes 
shifts the end-to-end distance statistics towards lower values 
of r. 



At this point, it is worth noting that £(£,#, u) is in- 
timately related to the moment generating function 
M(£) = (e« r ) = JdrG(r)e^ r . Using £(£,0,u) = 
£[S(r,0,u)] and[3l] 

G(r) = r 2 J d(cos8) J d 2 uG(r,6,u) (16) 

it is a straightforward matter to show that 

M(0 = dl [<*£($] , (17) 

where £(£) = J d(cos8) J d 2 u£(£, u). Hence, knowl- 
edge of £(£) directly provides us with the n th order mo- 
ments of the end-to-end distance r without explicit re- 
course to the RDF: 



(r») = d" +2 [e*£(0] 



(18) 



The RDF itself follows from £(£) by means of an inverse 
Laplace transform. Using eq. ( 16 1 and the definition of 
£(£) we find 



G(r) 



(19) 



where we used the linearity of the inverse Laplace trans- 



form. Inspection of eqs. (131 and (|14[) shows that £(£) 



contains poles of second order at {£p}peN- Thence we 
may invoke residues calculus in order to perform the in- 
verse Laplace transform (19 1 and write down G(r) as a 
sum over poles. Invoking our hitherto derived results we 
find to the desired order in the global modes 
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Here the following terms have been defined: 
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In eqs. (21) the sums 



S 



(fc) = 



E 



(21b) 
(21c) 

(22) 
(23) 



may actually be represented in terms of transcendental 
functions by means of an expansion in partial fractions 
(not shown). The function D(x), occurring in eq. (22) 
denotes Dawson's integral, D{x) = e~ x f Q dye v . 

As an aside we note that due to the fact that the expo- 
nential in eq. (21a) has a local maximum at q p = pn = c 
(c > 0), wavenumbers q p in the vicinity of the inverse 
Odijk length contribute the most to the sum in eq. 
(20). This is, of course, consistent with the observation 



that the polymer preferably stores length at wavenum- 
bers q R3 c = £7 . 

Repeating our calculations, stipulating rigid coin- 
cidence of polymer-fixed and channel-fixed reference 
frames — i.e. invoking hinged ends boundary conditions 
within E c — we are left with the following simplified ver- 
sion of the RDF 



Gh(r) = Z h 



OO 



p=l 



5p°\r), 



(24) 



which was determined before by Levi and Mecke in ref. 
|23j . Comparing eqs. ( [20| a nd (24), suggests to think of 
the RDF stated in eq. ( |20| as being made up of three 
structural parts: first, the coefficients Gp°\r), carrying all 
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FIG. 3: Mean end-to-end distance (r) as a function of con- 
finement strength c, calculated on the basis of G(r) (eq. ( |20[ ), 
solid line), Gh(r) (eq. (24 1, dashed line) and r 2 Qh(r) (dot- 
dashed line). The shaded area gives an illustrative measure 
of the impact of the energy optimizing effect. Inset: Only the 
first order corrected model exhibits a local minimum of (r) (c) 
at onsetting confinement. 



the statistical information related to the internal modes. 
Second, the common factor r 2 , and, third, the terms in 
curly brackets, both arising as a consequence of the global 
modes. Note that the global factor r 2 is simply a phase 
space factor and a consequence of the spherical integra- 
tions performed, whereas the non-trivial correction fac- 
tors in curly brackets are brought about by the explicit 
(#, u)-dependence of the confinement potential H po t- 

Depending on the confinement strength, the three 
structural parts have different impact on the chain's end- 
to-end distance statistics. Figure [2j comparing eqs. (20) 
and (24 1, illustrates this point. 



In the regime of vanishing (and weak) confinement, 
the global modes [32] are uniformly distributed, whence 
the non-trivial correction terms Gg) u p in eq. ( 20 1 vanish 



and G(r) reduces to the free-space solution determined 
in ref. |25j . In this context note that our expansion 
scheme reproduces the zero-confinement limit (where the 
expansion parameters become particularly large) even by 
zeroth order. Hence, our first order corrected model is 
supposed to give an excellent approximation for virtu- 
ally any choice of the confinement parameter c, safe for a 
very narrow regime of onsetting confinement, where the 
expansion parameters are still large and where the non- 
trivial effects of confinement are already perceptible. We 
will return to this point below. For now the important 
point to emphasize is that in the limit of vanishing con- 
finement the only fingerprint left behind by the global 
modes is embodied by the common phase space factor 
r 2 , which is lost by any approach where the positions 
and/or orientations of the polymer's ends are artificially 
constrained. From a physical point of view, the phase 
space factor r 2 arises for entropic reasons, since in 3 di- 
mensional embedding space the number of states consis- 
tent with a given contour {a„}, with end-to-end distance 
r, is proportional to the surface of a sphere with the 



same radius, and therefore to r 2 . Neglecting this factor 
enhances the probability of comparatively short configu- 
rations and thus shifts the mean end-to-end distance (r) 
below its actual value, which is exactly what can be ob- 
served in figs. [2ja) and [3] Apart from the simple factor 
r 2 , however, eqs. p0| ) and ( p4| are identical in the limit 
c — > 0, so that the shortcomings of Gh( r ) may easily be 
corrected. 

In the regime of intermediate confinement (fig. [2^t>)), 
eqs. (20) and (24) still deviate from each other, but this 



time the simple factor r alone is no longer capable of 
refining the behavior of Gh{r). Here, the essential cor- 
rections to be taken into account are contained in the 
nontrivial terms Gg) u p - In this case, the presence of the 
confining channel walls noticeably constrains the thermal 
undulations of the polymer. For any "frozen" contour, 
i.e. for any fixed set {a n }, however, the global modes 
allow for a shift of the whole contour with respect to the 
channel walls in such a way as to lower the potential en- 
ergy of the entire system. This energy-optimizing effect, 
brought about by the global modes and analytically con- 
tained within the nontrivial coefficients 0^9 , enhances 
the probability of relatively compact conformations. Ne- 
glecting these coefficients inevitably leads to an overesti- 
mation of the stretching effect due to the confining walls, 
which is precisely what can be observed in figs. [2jb) and 
[3] Moreover, note that since this effect mainly results 
from the topological restrictions imposed by the channel, 
the difference in chain statistics arising from the differ- 
ent conditions at the boundaries of the polymer is nearly 
independent of the chain's flexibility in this regime. 

Finally, once confinement is sufficiently strong, the 
global modes become so effectively constricted, that they 
finally lose their significance for the determination of the 
RDF, and the difference between G(r) and Gh( r ) van- 
ishes (not shown). Both functions approach the same 
limit, which is a delta-function at full extension. 

In addition to the quantitative differences discussed 
so far, the particular choice of boundary conditions also 
impacts the qualitative behavior of (r) as a function of 
the confinement strength. Figure [3] illustrates this point. 
Whereas hinging the polymer's ends to the channel axis 
implies a monotonically increasing end-to-end distance, 
accounting for the genuine conditions at the boundaries 
of the chain gives rise to a local minimum in (r)(c) at on- 
setting confinement. Although our first order corrected 
model is not capable of reproducing this feature in quan- 
titative detail, our approach unambiguously traces back 
this phenomenon — which has recently been reported in 
MC simulations [TH [20] — to the influence of the global 
modes. 

In conclusion, we gave an explicit, analytical account 
for the derivation of the radial distribution function 
(RDF), as well as for the n th order moments (r n ) of the 
end-to-end distance, for semiflexible polymers in cylin- 
drical confinements. Here we particularly focused on the 
use of proper boundary conditions and highlighted their 
influence within different confinement scenarios. To this 
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end, we introduced additional global modes in order to 
account for the possibility of the chain to perform global 
movements within the confining channel, and calculated 
their influence on the polymer's radial statistics to lead- 
ing order. In a regime of weak and vanishing confinement 
the global modes are trivially distributed and their effect 
on the RDF reduces to a simple phase space factor, which 
is of purely entropic origin. In contrast, in a regime of 
intermediate confinement the system's potential energy 
explicitly depends on the global modes, which therefore 
affect the RDF in a non-trivial way. We found that ne- 
glecting their influence in the latter regime inevitably 
results in an overestimation of the stretching effect due 



to confinement. Finally, in addition to these quantitative 
differences the global modes were found to be responsible 
for the formation of a local minimum in (r)(c) at onset- 
ting confinement. 
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